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Abstract 



We propose a regularization-independent method for studying a renormal- 
izable held theory nonperturbatively through its Dyson-Schwinger equations. 
Using QED4 as an example, we show how the coupled equations determining 
the nonperturbative fermion and photon propagators can be written entirely 
in terms of renormalized quantities, which renders the equations manifestly 
hnite in a regularization-independent manner. As an illustration of the tech- 
nique, we apply it to a study of the fermion propagator in quenched QED4 
with the Curtis-Pennington electron-photon vertex. At large momenta the 
mass function, and hence the anomalous mass dimension 7m,(a), is calculated 
analytically and we find excellent agreement with previous work. Finally, we 
show that for the CP vertex the perturbation expansion of 7m(tt) has a finite 
radius of convergence. 

I. INTRODUCTION 



As is well known, in studies of renormalizable relativistic field theories it is unavoidable 
to make use of a regulator in order to control the UV divergences of the theory, even though 
at the end this regulator is removed in order to obtain physical results. In perturbation 
theory the regulator and its removal via the renormalization procedure is usually little more 



* E-mail : akiziler ©physics . adelaide .edu.au 
''^E-mail: aschreib@physics.adelaide.edu.au 
E-mail : awilliam@physics . adelaide .edu.au 



1 



than an inconvenience, however in nonperturbative studies, which are largely numerical, 
much effort needs to made in order to perform this removal in a well-controlled manner. 

The best known example of the latter, of course, is lattice gauge theory [|I| where the 
regulator is the lattice spacing and its removal entails a careful numerical extrapolation to the 
continuum. Another popular tool, which is the subject of this investigation, is the numerical 
study of a truncated set of the Dyson-Schwinger equations of the theory where most 

often the regulator is simply a hard momentum cutoff A. Although these calculations are 
numerically far less demanding than the numerical integration of a functional integral, they 
are not entirely trivial either as they involve solving a set of integral equations for momenta 
p ranging over many orders of magnitude, from \p\ = to \p\ = A. Renormalization then 
requires a further numerical extrapolation, with fixed boundary conditions provided by the 
relevant renormalization scheme, to A — oo. 

That this removal of the regulator can be done numerically was demonstrated explicitly 
in Refs. [0-0] within subtractively renormalized quenched QED, using an off-shell renor- 
malization scheme. However, as emphasized in those works and pointed out first by Dong, 
Munczek and Roberts ^ (it was also apparently realized at an earlier stage by Penning- 
ton 0), the use of a momentum cut-off as a regulator is not entirely satisfactory as it violates 
translational invariance (in momentum space) which leads to a violation of gauge invariance 
in the final results, even after A is taken to infinity. More precisely, this violation of gauge 
invariance was observed in quenched QED calculations employing the Curtis-Pennington 
(CP) electron-photon vertex and may be traced back to a certain, logarithmically diver- 
gent, 4-dimensional momentum integral which vanishes because of rotational symmetry at 
all A < oo, but leads to a finite contribution for A — > oo. It is this discontinuous behaviour 
as a function of A which complicates correct numerical renormalization with this regulator. 
Incorrect results will be obtained unless care is taken to drop 'gauge covariance violating 
terms'. Although this recipe appears to work in practise (see the Appendix of Ref. as 
well as Ref. for a more detailed discussion) it is fair to say that it is ad hoc and it would 
be more desirable to avoid it altogether. 

Dimensional regularization, of course, does not break translational invariance and hence 
will not lead to the spurious gauge covariance violations outlined above. However, as shown 



in Ref. |TT| , the use of this regulator in nonperturbative studies of chiral symmetry breaking 



introduces a new unwelcome feature: In a dimensionally regulated theory a fermion mass 
will be generated dynamically at all values of the coupling constant a in D ^ 4 dimensions 
for any vertex which only does so above a finite critical coupling a^r in D = 4 dimensions. 
In other words, again discontinuous behaviour (this time as a function of A — D) prevents 
a naive removal of the regulator. Ultimately, it is this feature which limited the accuracy 
with which one could determine the critical coupling associated with the CP vertex within 
dimensional regularization in Ref. Jill. ■'^^ addition, of course, with this regulator the inte- 
gration range now extends to \p\ —>■ oo, compounding the numerical difficulties mentioned 
at the start. In short, although it was shown in Ref. that the results obtained within 
dimensional regularization agree with those obtained by employing a cut-off within some 
limited precision, the numerical difficulties encountered in implementing this regulator has 
limited its further use. 

It is the purpose of this letter to point out that all the difficulties mentioned above may 
be circumvented by removing the regulator analytically rather than numerically. This may 



2 



be done by formulating the Dyson-Schwinger equations in terms of renormalized quantities 
only, where the dependence on the mass scale introduced by the regulator is traded for the 
momentum scale /i at which the theory is renormalized before performing any numerical 
calculations. This way dimensional regularization (or any other regulator which doesn't 
violate gauge covariance) can be used and, more importantly from a numerical point of view, 
removal of the regulator means that no longer does one have to solve integral equations 
involving mass scales of vastly different orders of magnitude (and then, in addition, take 
a limit in which one of the scales goes to infinity). Rather, the important scales in the 
problem become scales of 'physical' importance, such as the renormalized mass and 
the renormalization scale fi at which this mass is defined. It is therefore to be expected 
that the dominant contributions to any integrands will be from a relatively small region of 
momenta. This feature is generic; i.e. it is independent of the particular vertex Ansatz that 
one makes use of and it remains a valid consideration for an arbitrary renormalizable field 
theory. Hence, although we shall illustrate this approach only within QED (and numerically 
only within quenched QED), it is our hope that ultimately this 'regularization- free' approach 
will prove its utility elsewhere. 

To a certain extent there is a price to be paid for using the approach advocated above: one 
looses all contact to the bare theory. In particular, we stress that because of this one cannot 
study dynamical chiral symmetry breaking in this approach by simply setting the bare mass 
mo equal to zero and investigating at what point a dynamical fermion mass is generated. It 
is in this crucial point that we differ from the analysis in Ref. [|1^] where, within quenched 
QED, a removal of the regulator was attempted in a manner which has some similarity 
to what we do here. Indeed, as emphasized particularly by Miranksy (Refs. [1^, as 
well as Ch 10.7 of Ref. 0), some care needs to be taken with the treatment of the bare 
mass while removing the regulator in order to avoid drawing incorrect conclusions about 
the presence or absence of dynamical chiral symmetry breaking in gauge theories. In the 
present approach this problem is avoided: chiral symmetry breaking, both explicit as well 
as dynamical, is characterized by a nonzero renormalized mass M^, and hence this quantity 
in itself does not distinguish between these two possibilities. In Section III, we shall take 
(within quenched QED) the somewhat more indirect signal of the appearance of oscillations 
in the renormalized fermion mass function [|l^] as an indicator of the onset of dynamical 
chiral symmetry breaking. 

This paper is organized as follows: in the next section all reference to the bare Lagrangian 
will be explicitly eliminated by re-expressing the Dyson-Schwinger equations of full QED in 
terms of renormalized quantities only. As a first application of this seemingly innocuous step 
we then re-examine, in the second part of this paper, the large momentum behaviour of the 
electron propagator in quenched QED4 (employing the Curtis-Pennington vertex [jl2[)- In 
particular, because of the absence of a cut-off, scaling invariance in this momentum region 
is restored in the equation for the electron's mass function and hence this equation can be 
solved analytically. We compare the results with existing numerical studies and also use 
it to derive the analytic form of the anomalous mass dimension 7m(tt) within this offshell 
momentum subtraction scheme. Finally, we discuss the analytic behaviour of 7m(«) as a 
function of the coupling and conclude. 
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II. REGULARIZATION-INDEPENDENT FORMULATION 



The Dyson-Schwinger equations for the renormahzed electron and photon propagators 
(Fig. |1|) are well known. For example, employing dimensional regularization, we have 

S-\fi;p) = ZMS'-\p) - zZ,{f,)e,'J r"(/x; p, fc) ^(/x; fc)7^ g) (2.1) 

^ZMS'-\p) - Zi(/i)S(p) 

D~}{li]q) = Z^{^i)Dlp (q) + iZi{fi)e^^NF tr J r„(/i; p, fc) 5(/i; fc) 7^ 5(/i; (2.2) 

^ Z,ifi)Dl^-\q) + Zi(^)n„^(g) 
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FIG. 1. The Dyson-Schwinger equations for the fermion and photon propagators in QED 

Here Np is the number of species of fermions in the theory and the constants Zi{fi), Z2{fi) 
and Z^lfi) are the vertex, fermion wave function and photon wave function renormalization 
constants, respectively. We have assumed, for convenience, that the Nf fermions are identical 
(i.e. they have equal mass), the generalization being straightforward. S{^;p) and Da/s^fi; q) 
are the full fermion and photon propagators renormahzed at the momentum scale /i {S^{p) 
and D^p{q) are their bare counterparts), while the full (proper) fermion-photon vertex is 
rQ,(/i; k,p). These propagators are usually expressed in terms of scalar functions Zi^ji^^p^), 
Mip^) and (^(/x^; q^) through 



S{fi;p) 



Z{fi';p') 

i) - M{p^ 



q2 



GifM ■,q ){9al3 



g2 



(2.3) 



where ^ is the renormahzed covariant gauge fixing parameter, so that the fermion's equation 
separates into two scalar equations 



Z-\n';p') 
M{p')Z-\^'-p') 



Z2(yU) - Zi{n) Sd(p^) 

Z2(/i) mo + Zi{n) T.sip'^) 



(2.4) 
(2.5) 
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Here and Ss(p^) refer, respectively, to the Dirac odd and even parts of i.e. 

=i) ^d{p'^) + Ss(p^). The bar over these quantities indicates that we have exphcitly 
separated out the renormahzation constants Zi{fi); also, note that we do not indicate the 
imphcit dependence of l^d^sip^), through the functions Z{ii^-,p^) and Daf3{fi;q), on /x^. The 
equivalent equation for the photon wavefunction renormahzation g^) is given by 

G-\f^';q') = Z.ifi) + Z,if^)n{q') , (2.6) 

where n"^(g) = — (^g°'^ — n(g^). Note that, as we are free to use a regulator which 
respects gauge covariance and will use a vertex which satisfies the Ward-Takahashi identity, 
n°^(g) does not contain spurious quadratic divergences |]15| and is transverse as it should be. 



Hence the term involving the gauge parameter ^ on the left hand side of Eq. ( p.2| ) precisely 
cancels the corresponding term involving in D^jj on the right hand side (i.e. = ^sO- 
The renormahzation constants Zi{fi) as well as mo are fixed by the boundary conditions 

Z(/i2;^2) _ ^ ^(^2) ^ Gifi^;fi^) = 1 . (2.7) 

(Although we shall not do so here, in practice the fermion and photon can be renormalized 
at two different scales /i^ and /i'^. In QED it is usual to choose the latter as /i'^ = 0, in which 
case e^'^/{4:7i) just becomes the physically measured fine structure constant a.) Because of 
the Ward-Takahashi identity, which should be satisfied by the vertex Ansatz for Fq,, one has 

Zi(/i) = Z2{fi). 

In order to avoid cumbersome notation, we have not explicitly indicated functional de- 
pendence on the regulator in Eqs. ( p.l| )-(p^). Most of the above quantities, i.e. Z'(/i^;p^), 
MijP') and G(/i^; g^), as well as Zi_2,3(/u), are functions of the regulator. Indeed, if one keeps 
and fixed as one removes the regulator, which is the standard textbook prescription, 
the bare fermion mass mg and the bare charge also become regulator dependent. 

As one removes the regulator, the integrals on the right hand side of Eqs. ( |2.1|) and ( p. 21) 



and hence Erf(p^), Ss(p^) and n(g^) diverge logarithmically (or, to be more precise, in 
dimensional regularization they develop singularities at D = 4). It is the defining feature 
of a renormalizable field theory that these divergences may be absorbed into the constants 
^i,2.3(yu) and into the bare mass mo, rendering finite limits for Z{^'^\ p^), M(p^) and G(/i^; q^). 
However, we can make use of the independence of Zi 2,3(Ai) and mo in order to eliminate 
these constants from the above equations. For example, by evaluating Eq. ( |2.4| ) at a second 
momentum which, for convenience, we take to be = /i^, and forming an appropriate 
difference one obtains 

Zif^';p') = 1 + Zif^';p')i:,ip') - S,(/i2) (2.8) 



where we have made use of the renormahzation condition (|2.7| ). Because the left hand side 
of this equation is finite as the regulator is removed, the right hand side must also be (even 
though, of course, the individual terms on the RHS separately diverge). In a similar way 
one obtains 

M(/) = + [M(p2)S,(p2) + S,(p2)] - [M^Ildip^) + S,(/x2)] (2.9) 

and 
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G 



1 + 



G 



(2.10) 



We note in passing that the perturbative solution of these three equations defines 'finite 

^ We 



QED' in a rather different sense to that developed by Johnson, Baker and Willey |16 



do not pursue this further here; rather, the main point of this letter is to point out that 
these renormalized equations ( |2.8| ) to ( p.lO| ) (with the regulator removed, of course), or 
equivalently their counterparts in QCD or any other renormalizable field theory, provide a 
starting point for nonperturbative investigations which has some significant advantages over 
the usual treatment found in the literature. To illustrate this we now turn to the well-studied 
example of quenched QED. 



III. APPLICATION TO QUENCHED QED4 



For quenched QED Eq. (|2.10|) becomes irrelevant and Eqs. (|2.8|) and ( p.9|) reduce to 
one-dimensional integral equations. In particular, if we employ the CP vertex (which, 
importantly for the present approach, satisfies the Ward-Takahashi identity and respects 
renormalizability) and make use of a regulator which does not violate gauge covariance then 
these renormahzed equations become, in Euclidean space. 



4tt 



+ 



a 



[A;2 + M2(P)] 
°° dp 



4tt Jo [A;2 + M2(A;2)] 
r/A-2 



a 

^^4^70 [A;2 + M2(A;2)] 



Z(/i2;fc2) 

Z(/i2;p2)/(A:^/) - lik^fj.'] 
J{k\^^) 



+ M{p')I{k\p')-M^I{k\^i^ 



(3.1) 



(3.2) 



where I{k ,p ) and J{k ,p ) are those of Table. 1 of Atkinson et al.'s paper |]T^, apart from 



the term proportional to ^ (^1 

+ £i^J in J(fc2,p2) in their paper which originated from the 

violation of gauge covariance by their cutoff mentioned in the introduction. As promised, 
these equations are finite in the absence of a regulator. The equation for the wavefunction 
renormalization function (Eq. |3.1| ) is the same as Eq. (10) in Ref. |jl3|, however our treatment 
of the mass function differs from theirs. We shall not solve these equations numerically here 
but, rather, shall make use of the fact that, for momenta k"^ larger than M'^{k'^), the absence 
of a cut-off implies that they become scaling invariant. Hence we can solve them analytically 
in this region for finite M^, in the same way as was done by Atkinson et al. fl^ for the 
chirally symmetric theory (i.e. Mq = = 0). 



A. Large Momentum Behaviour in Quenched QED4 



For A;2 » M2(A;2) 
linear terms in M{k'^) 



Eqs. (|3.1j , pTZj ) may be expanded in powers of 
we obtain 



fc2 • 



Keeping at most 
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ai f^'' dk'^ 2 r 



(3.3) 



2^ < f dk^zi^i^-e) ,,,,2> 



47r [jo Z{ij?]p'^ 



f^'dl^ 

/i^ 



Z{i^';k')M{k') 



3a /■°° ,,2 
Stt io 



2M{k 



2 g(A'";fc^) 



Py p"^ — k"^ 



Zifi^; e) p\ M{p^) - M{k^ 
Z{^'^;p'^)p'i, p"^ — k"^ 




(3.4) 



where p\ (p^) is the minimum (maximum) of p"^ and k'^. 

Because of scahng invariance, these equations are solved by a power law Ansatz for 
Z{fi^;p'^) and M{p^), i.e. 



Z{f^';p'' 



p_ 



M{p^ 




(3.5) 



The reader should note that the overall scale of M{p'^) is fixed by the renormalization 
condition. Explicitly, one finds u and s to be given by 



47r 



(3.6) 



3u 
2^ 



2tt cot STT — IT cot z/vr + 37r cot(z/ 



S)TC + 

3 



2 11 



1 



[ly - s) 



Z/-S + 1 



+ 



V — s 



(3.7) 



These are the same powers as those found in Ref. [[T^] (corrected for the gauge covariance 



violating term) . We emphasize that in contrast to the unrenormalized equations in Ref. [jT7 



for the renormalized equations this power Ansatz is valid whether or not the bare mass mo is 
zero. We demonstrate this explicitly in Fig. (^, where the above asymptotic form of the mass 
function (dashed line) is compared to previously published numerical work (Ref. 0) with 
mo 7^ (solid line). As one can see there is excellent agreement in the region p^ » M^(p^). 
Quantitatively, one sees in Table [1] that the logarithmic slope s calculated from Eq. ( p.7| ) for 
a variety of gauge parameters and couplings agrees, in most cases, with the same quantity 
extracted from the published numerical work to better than 1 part in 10^. 

We note that in the above calculations a < Ucr- There are, in general, several real 
solutions s to Eq. ( p.7| ) at these values of the coupling (see also Ref. [0). However, only 
one of these matches smoothly onto the perturbative solution and it is this one which 
is used above. The other solution appears to be a spurious byproduct of the linearized 
approximation (|3.4|) to Eq. (|3.2|) (Solutions which do not match smoothly onto perturbation 
theory can arise in Eq. ( |3.4| ) if the integrals diverge like - due to IR divergences; this cannot 
occur in the exact Eq. ( p.2| ) because of the regulating mass in the denominator). Because s 
is real there are no oscillations in the mass function and hence there is no chiral symmetry 
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breaking for these couplings. For a > acr, however, the solutions s become complex. As 
Eq. (|3.7| ) defines s as a real analytic function of a, it is clear that if s is a solution, then so is 
s*. In this case the correct asymptotic tail to the oscillating, but real, mass function (Eq. p.2| ) 
is obtained by a suitable superposition of these two particular solutions to Eq. ( p.4|) . The 
oscillations are the indication that, above acr, chiral symmetry is broken in quenched QED 
employing the CP vertex. 
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TABLE. I. Comparison of asymptotic behaviour 
of mass function (Eq. (3.5)) between the numerical 
work [5,6] and Eq. (3.7) for a given a and 5. 





Q 


.s 

Xumerical[5] 


Eq. '(3.7) 





0.1 


0.0244890 


0.0244888 





0.2 


0.0504181 


0.0504176 





0.3 


0.0781811 


0.0781793 





0.4 


0.10832.56 


0.1083229 





0.5 


0.1416687 


0.1416653 


0.25 


0.5 


0.143026 


0.1430124 


0.5 


0.5 


0.144349 


0.1443495 





0.6 


0.1795307 


0.1795287 





0.7 


0.2243289 


0.2243399 





0.8 


0.2816516 


0.2816491 


0.25 


0.8 


0.2874036 


0.2875670 


0.5 


0.8 


0.2933384 


0.2935501 



FIG. 2. Comparison of a numerical solution for the mass function M(p^) (solid line, obtained 
from Ref. Q) and the analytical expectation for its asymptotic behaviour (dashed line, Eq. 3.7). 
The renormalized mass is = 400 (arbitrary units), while a = 0.6 and ^ = 0.5. 



B. The anomalous mass dimension 

The analytical expression for the asymptotic behaviour of the mass function is closely 
related to the anomalous mass dimension of the theory, i.e. 

dlogZm „ 

7m = —p. ^ 2s , (3.8) 

dlog/i 

where we have used = tuq / and the asymptotic form of the mass function in Eq. 



Hence Eq. ( |3.7|) defines, within the off-shell renormalization scheme used here, the prediction 
for the asymptotic value of the anomalous mass dimension obtained with the CP vertex as a 
function of the couphng and the gauge parameter ^. For convenience, let us restrict ourselves 
to Landau gauge, in which case one obtains 
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It is amusing to note the similarity of this equation, obtained within the framework 
of Dyson-Schwinger equation studies, to that obtained for the anomalous mass dimension 
within the MS scheme using a completely different nonperturbative approach JlSl, namely 



one based on the use of a variational principle for the evaluation of the functional integral, 
i.e. 

« = ^ (1 + cot f ^ . (3.10) 



3 ^ ' V2 1 + 7, 

Not only are both equations implicit equations for 7^, but even the same transcendental 
function is involved in both. Actually, the similarity does not end there. It is a simple matter 
to numerically solve Eq. (|3.9|) for complex a in order to investigate its analytic structure. 
The result, for both the imaginary and real part of the anomalous mass dimension, is shown 
in Fig. d^). A branch point in a is clearly visible on the positive real axis (the direction of 
the cut is arbitrary). Its position may be determined via a bifurcation analysis of Eq. ( p.7[ ) 
or, equivalently, Eq. ( p.9|) . Its numerical value is a = 0.933667 (in Feynman gauge it is 
a = 0.890712) and it coincides, of course, with the well-known critical coupling associated 
with the CP vertex |lT7| , p!T| above which the theory breaks chiral symmetry. The variational 



result for the same quantity in the MS scheme, Eq. (|3.10|) , also results in a branch point, 
however this time it is at a complex value of a |T8| , p!9| . Hence no sign of chiral symmetry 
breaking for a > was seen in 7™'' in that work. 
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2 Im 7^ 
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-0.2 




Im a 



FIG. 3. The real and imaginary part of the anomalous mass dimension versus complex coupling 
constant. 



Apart from the presence of a branch point at finite a, the other notable feature in 
Fig. (^ (also found for 7™ (a) in Refs. [|^,|^) is the absence of a branch point, or any other 
nontrivial analytic structure, at a = 0, i.e. the radius of convergence of the perturbation 
expansion of 7m(Q;) as obtained with the CP vertex is finite. This is in marked contrast to 



general expectations from studies of high order perturbation theory |2^, where the factorial 
growth in the number of diagrams results in a vanishing of the radius of convergence of the 
perturbation expansion of a general Green function. It is difficult to specify whether or not 
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the CP vertex correctly reproduces the factorial growth in the number of diagrams. What 
can be said with certainty, however, is that clearly the Dyson- Schwinger equation for the 
fermion propagator does not itself generate a branch point at a = like it does, for example, 
at a = Ocr- If a branch point at a = is present in the exact theory then it already needs 
to be present in the vertex itself. 

Finally, associated with the above analj^ical structure, one can derive the behaviour of 
high order perturbation theory for 7™. (a) from Eq. (^.91), using the methods discussed in [jiy 
We define the expansion coefficients c„ through 



7m (a 



) = E 



n=l 



a 

a = — 

TT 



(3.11) 



so that the radius of convergence (in a) is given by lim„^oo |cn|^^^"- The asymptotic be- 
haviour of these coefficients may be derived most easily by converting Eq. ( p.7|) into a 
differential equation: 



3a2 



2^ + 



-— + l + 3/(a)+j(o.) 



7m(a) + 3/'(a) +^'(a) 



(3.12) 



where we have defined /(a) = 2/7^(0) and g{a) = — 7m(a)/2]. These functions have 
the perturbative expansions 



r) 00 

— E^n a" 



J2 a" 

n=0 



where the coefficients 6„ and dn are related to the c„'s through PT 



(3.13) 



bo = 1 



do 
dn. 



k=l 



Cl 



1 

1 " 
^ fe=l 



n > 



(3.14) 



These equations, together with Eqs. ( p.ll| ) and ( |3.13|) substituted into Eq. ( |3.12| ), define the 
expansion coefficients Cn and may be solved iteratively to very high orders. For the first few 
orders we obtain 



7m(a) 



3 a; 9 /a 

2^ ^ 8 W 
81 



a 



TT 



256 



(e + TT^ 



2560 



'105 + 207r^ 



TT 



while for large n the behaviour of the coefficients approaches the form 



a 



cr 
TT 



-n g— /3o 



3/2 



(3.15) 



(3.16) 



with /^o ~ 0.52. The absence of a term of the form n" indicates the absence of factorial 
growth mentioned previously and the power 3/2 in the denominator is indicative of the fact 
that close to the branch point the cut has the character of a square root. 
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IV. CONCLUSIONS AND OUTLOOK 



In this letter we have argued that there are significant advantages in performing nonper- 
turbative calculations based on the use of Dyson-Schwinger equations in a manner different 
to that commonly used in the literature. The central idea is to work with renormalized 
quantities only, allowing one to completely eliminate the regulator. Within QED, this al- 
lows one to avoid spurious problems associated with the usual cut-off regulator, such as the 
breaking of gauge covariance in the fermion propagator and the introduction of unwanted 
quadratic divergences into the photon propagator. Furthermore we have shown that, within 
quenched QED, it allows one to gain an analytical understanding of both the mass function 
at large momenta as well as the analytic structure of the anomalous mass dimension of this 
theory. We expect that the main advantage of the method, however, will be that one can 
restrict the numerical task of solving the Dyson-Schwinger equations of any renormalizable 
theory to a relatively small region of momenta. For example, we anticipate that this feature 
will make it feasible to study full, rather than quenched, QED4. Furthermore, it is our hope 
that the numerical task of solving a coupled set of Dyson-Schwinger equations is sufficiently 
simplified by the approach discussed here to allow one to consider pushing the truncation 
of this set beyond the equations for the two-point functions. 



ACKNOWLEDGMENTS 

We are grateful to F. Hawes for providing us with the data of Ref. used in preparing 
Table 1 and T. Sizer for helpful conversations about this work, as well as a careful reading 
of this manuscript. 



11 



REFERENCES 



[1] H. J. Rothe, Lattice Gauge Theories: an Introduction, (World Scientific, Singapore, 
1992). 

[2] P. I. Fomin, V. P. Gusynin, V. A. Miransky and Yu. A. Sitenko, Riv. Nuovo Cim. 6 
(1983) 1. 

[3] V. A. Miransky, Dynamical Symmetry Breaking in Quantum Field Theories, (World 

Scientific, Singapore, 1993). 
[4] CD. Roberts and A.G. Williams, Prog. Part. Nuc. Phys. 33 (1994) 477. 
[5] F.T. Hawes and A.G. Williams, Phys. Rev. D 51 (1995) 3081. 
[6] F.T. Hawes, A.G. Williams and CD. Roberts, Phys. Rev. D 54 (1996) 5361. 
[7] F.T. Hawes, T. Sizer and A.G. Williams, Phys. Rev. D 55 (1997) 3866. 
[8] Z. Dong, H. Munczek, and C D. Roberts, Phys. Lett. 333B (1994) 536. 
[9] See Footnote 1 in A. Bashir and M. R. Pennington, Phys. Rev. D 50 (1994) 7679. 
[10] A. Kizilersii, T. Sizer and A. G. Williams, ADP-99-43-T380, |hep-ph/000lTT7| . 
[11] A.W. Schreiber, T. Sizer and A.G. Wilhams, Phys. Rev. D 58 (1998) 125014 ; V. P. 

Gusynin, A. W. Schreiber, T. Sizer and A.G. Wilhams, Phys. Rev. D 60 (1999) 065007. 
[12] D.C Curtis and M.R. Pennington, Phys. Rev. D 42 (1990) 4165 and references therein. 
[13] D.C Curtis and M.R. Pennington, Phys. Rev. D 46 (1992) 2663. 
[14] V. A. Miransky, Phys. Lett. 165B (1985) 401; ibid, II Nuovo Cim. 90A (1985) 149; see 

also D. Atkinson and P. W. Johnson, Phys. Rev. D35 (1987) 1943. 
[15] K.-I. Kondo, H. Mino and H. Nakatani, Mod. Phys. Lett. A17 (1992) 1509; 

J. C. R. Bloch and M.R. Pennington, Mod. Phys. Lett. AlO (1995) 1225. 
[16] K. Johnson, M. Baker and R. Willey, Phys. Rev. 136 (1964) Bllll. 
[17] D. Atkinson et al., Phys. Lett. 329B (1994) 117. 

[18] C. Alexandrou, R. Rosenfelder and A.W. Schreiber, Phys. Rev. D 62 (2000) 085009. 

[19] A.W. Schreiber, R. Rosenfelder and C. Alexandrou, in "Proceedings of the 3rd Inter- 
national Symposium on Symmetries in Subatomic Physics", eds. X.-H. Guo and A. 
W. Thomas, to be published by AIP Conference Proceedings, ADP-00-20 /T403, [hep- 

I th/000718^ . 

[20] L. N. Lipatov, Sov. Phys. JETP 45 (1977) 216; for a review, see E. Bogomol'nyi, V. 

A. Fateev and L. N. Lipatov, Physics Reviews 2 (1980) 247 as well as the collection of 

papers in "Large-order behaviour of perturbation theory", eds. J.C Le Guillou and J. 

Zinn- Justin, North-Holland (1990). 
[21] I. S. Gradshteyn and I. M. Ryzhik: Table of integrals, series and products, 4th edition, 

Academic Press (1980), p. 14. 



12 



